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Abstract 



The electronic structure of the helium atom in the magnetic field regime B — 0— 100a.it. is investigated, 
oo; _ ... . . ... . 

^\ ' using a full configuration interaction approach which is based on a nonlinearly optimized anisotropic 

o ' 

for the magnetic quantum number M = — 1 and for both even and odd z-parity as well as singlet and 



Gaussian basis set of one-particle functions. The corresponding generalized eigenvalue problem is solved 



triplet spin symmetry. Accurate total electronic energies of the ground state and the first four excitations 
in each subspace as well as their one-electron ionization energies are presented as a function of the magnetic 
Oh' field. Additionally we present energies for electromagnetic transitions within the M = — 1 subspace and 

between the M = — 1 subspace and the M = subspace treated in a previous work. A complete table of 
> ' wavelengths and field strengths for the detected stationary points is given. 



1 Introduction 

Motivated by the astrophysical discovery of strong magnetic fields on the surfaces of white dwarfs (< 10 5 
Tesla) and neutron stars (~ 10 8 Tesla), the behaviour and the properties of matter in strong magnetic 
fields has increasingly attracted interest. The theoretical description of atoms in strong magnetic fields 
is well covered in the literature only for the case of the hydrogen atom (see refs.[[j]-|7| and refs. therein). 
Until recently our knowledge about atoms with more than one electron in a strong magnetic field has been 
relatively sparse and definitely not sufficient for a comparison of the corresponding theoretical data with the 
mysterious absorption edges f8|HlO| in the spectrum of the magnetic white dwarf GD229. Those have for the 
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first time been identified as helium lines by Jordan et a^Ol m 1998. This overwhelming evidence was based 
on results of highly accurate ah initio calculations performed by Becken and Schmelcher in 1998, a part of 



which has been published in ref . 12]- For a more detailed overview over the various theoretical approaches 
to the helium atom in a strong magnetic field and the corresponding literature before 1998 we refer the 



reader to [12| and in particular the references therein. In (T^| a fully correlated configuration interaction 
approach has been applied to the helium atom in a strong magnetic field. Total energies for spin singlet 
and triplet states for both positive and negative z-parity in the subspace of vanishing magnetic quantum 
number M = have been provided, thereby covering the regime of magnetic fields strengths from B = to 
B = 100a.it. [B = la.u. corresponds to 2.35 • 10 5 Tesla). Additionally all the transition energies within the 
M = subspace have been presented and discussed there, including the stationary components with respect 



to the field dependence which have been the key tool for the comparison with the observed spectra [11]. 

The aim of the present paper is to provide important data for states with the magnetic quantum number 
M = — 1. Some data for states with positive z-parity and M = — 1 have already successfully been used for 



a comparison with the astronomical observation(in [11], we used the corresponding stationary transitions 
in a certain field and wavelength regime). The complete data are presented here for the first time. Data 
resulting from states with negative z-parity are very recent results. The latter further extend our knowledge 
on the helium atom in a strong magnetic field and permit to investigate additional transitions. 

For the case of triplet spin symmetry, Jones et al |13| very recently used a released-phase quantum Monte 
Carlo method for calculating accurate data. However, they cover only three field strengths, investigate less 
excited states for common symmetries and do not study the spin singlet states at all. Nevertheless for the 
common values of the field strength they confirm our data to several digits. In contrast to this the energies 
presented by Scrinzi fl4}| , though in principle variational, are in the presence of a magnetic field significantly 
lower than ours. Carefully comparing the energies, it appears to us that they are systematically too low, 
i.e. most probably due to a numerical error. 

We remark that performing calculations for finite magnetic quantum numbers requires within our ap- 
proach drastically improved computational techniques (in comparison to the case M = 0) for keeping the 
CPU time affordable. A summary of these techniques is provided in the appendix of the present paper. 

The starting point of the present paper is the nonrelativistic Hamiltonian of the helium atom with 
infinite nuclear mass in a magnetic field as given in section 2. To be self-contained we briefly discuss 
the Hamiltonian's symmetries and provide a description of the basis set as well as the full configuration 
interaction approach (for more details see ref.[Q). We introduce a maximal set of conserved quantities, 
chosen to be the total spin S 2 , the z-component S z of the total spin, the total spatial magnetic quantum 
number M and the total spatial z-parity Tl z . These symmetries serve for classifying the results for the 
energies for M = —1 in sec. 3. In each of the two subspaces for positive and negative z-parity we present the 
total energies and the ionization energies of the ground state and the first four excited states for singlet and 
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triplet spin symmetry. Additionally we consider in sec. 3 all the transitions within the M = — 1 subspace as 
well as all the possible transitions between the M = — 1 states treated in the present paper and the M = 
states given in ref.[Q. The wavelengths of all the stationary components are provided, being the basic 
ingredient for the successful comparison of theoretical data with the spectrum of magnetic white dwarfs in 
general and specifically for GD229. 

2 Hamiltonian, Symmetries and basis sets 
2.1 Hamiltonian and Symmetries 

Assuming the magnetic field to point in the +z-direction, the Hamiltonian reads 

1 



H = ^ P l + lBl zi + Bl( X * + yt)-^ + Bs zi ) + ^—-^ (1) 

The one-particle operators in eq.(0) are the Coulomb potential energies — JL of the electrons in the field 

of the nucleus as well as their kinetic energies, here splitted into the parts ip?, the Zeeman terms ^Bl zi , 

the diamagnetic terms ^-( x i + Vi ) an d their spin energies Bs Z i. The electron-electron repulsion energy is 
represented by the two-particle operator | r2 ^ ri | ■ We remark that we use an electron spin ^-factor equal 2, and 
any more accurate value for it can be simply incorporated by shifting the final total energies correspondingly. 
For remarks on the influence of relativistic effects and on a scaling relation taking into account the finite 
nuclear mass we refer the reader to ref.[||, 12, D|. 



Analogously to ref.[12|, we exploit that there exist four independent commuting conserved quantities: 
the total spin S 2 , the z-component S z of the total spin, the z- component L z of the total angular momentum 
and the total spatial z-parity Ii z . 



2.2 Basis sets 

For constructing a two-particle basis set of eigenfunctions of the above mentioned conserved quantities, our 
central ingredient is an anisotropic Gaussian basis set of one-particle functions 

$i(p,<p,z) = p n p, z n ^ e - aip2 -^ z2 e im ^ i = l,... ,n , (2) 

which are themselves eigenfunctions of the corresponding one-particle operators of the mentioned conserved 
quantities. The parameters n pi and n zi are restricted by 

n pi = \rrii\ + 2k{ ; ki = 0, 1, 2, ... with m, = ... — 2, — 1, 0, 1, 2, ... (3) 
n zi = ir zi +2U; k = 0,1,2,... with n zi = 0, 1 (4) 
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whereas the nonlinear variational parameters a% and /3j are positive and have to be nonlinearly optimized 
for each field strength as described in ref.[12|. For each one-particle subspace of given symmetry we used 



an algorithm for determining the nonlinear parameters aj and [3i such that the states of the hydrogen atom 
or the He + ion for that symmetry were optimally described. We emphasize that this procedure gives rise to 
considerable effort since it has to be repeated for each field strength separately. 
We construct a basis set of spatial two-particle states by 

\ip q ) := b\b)\G) i = l,...,n , j = i,...,n , (5) 

where b\ is the creation operator of the i-th one-particle state \i) = b\ |0) whose position representation is 
given by eq.(||). The spin space is spanned by spin singlet or spin triplet states, and therefore the operators 
b\ have to be chosen bosonic or fermionic, respectively. Selecting combinations with 

mi + mj = M , mod(ir zi + ir Z j, 2) = IT^ , (6) 

we achieve the two-particle states (|5|) to be a basis set within the subspace for given total symmetries M 
and Tl z . The number N of two-particle basis states is thus in general smaller than n(n + l)/2. 

We perform a full Configuration Interaction (full CI) approach by representing the Hamiltonian in 
a basis whose spatial part is given by the in general nonorthonormal states (||). Since the spin part BJ2 s zi 
of the Hamiltonian can trivially be taken into account by a shift of the energies it is sufficient to represent 
the spatial part of the Hamiltonian H and the overlap S by 

Spq = (lpp\lpq) , H pg = (lp p \H\lp q ) (7) 

The matrices S and H are Hermitian, and the overlap S is additionally positive definite. Furthermore the 
matrix elements turn out to be real. The finite-dimensional generalized real-symmetric eigenvalue problem 

(K-E§)-c = (8) 

provides eigenvalues E which are variational upper bounds to the exact eigenvalues of the Hamiltionian (|l]) 
within each subspace of given M and H z . 

2.3 Matrix elements 

For calculating the matrix elements of the spatial part of the Hamiltonian (|]), we rewrite the former in 
second quantization, H = Hj + Hu, where Hj and Hu denote the second-quantized counterparts of the 
familiar one- and two-particle operators whose position representations read 

F 7 ( P ,r) = ipS+iB-l + pV + y 2 )-^ iT//(r!,r 2 ) = (9) 
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Now, with \tp q ) := b\bj |0) and \ip p ) := b\b\ |0) a straightforward calculation leads to 

{^ p \%) = (i\k) ± <t|i> (j|fc> (10) 

= <i|S-j|Ar> 010 ± (i|H 7 |Z) (j\k) 

+ (j\Hj\l) (i\k) ± (j\Hj\k) (11) 

il> p \Hu\iJ> q ) = (ij\H n \kl) ± (ij\H n \lk) , (12) 



where \ij) := \i) ® |j) and where the sign '±' stands for '+' in the singlet case and for '— ' in the triplet 
case. 

For the relatively simple evaluation of the n(n + l)/2 different one-particle overlaps (i\k) and matrix 



elements (i\Hj\k) we refer the reader to appendices A,B in ref.[12|. The two-particle matrix elements 
(ij\Hjj\kl) are by no means trivial, in particular in view of the fact that their accurate and fast evaluation is 
necessary in order to build up the Hamiltonian matrix in an affordable amount of CPU time. In ref.|l2| we 
discussed a method using a decomposition in Cartesian coordinates which expresses the two-particle matrix 
elements in series of hyper geometric functions whose evaluation has been performed by highly efficient 
analytical continuation formulas. The latter are necessary in order to keep the CPU time acceptable since 
the number of different two-particle matrix elements is of the order N(N + l)/2 rather than n(n + l)/2. 
However, the Cartesian decomposition becomes more and more inefficient with increasing magnetic quantum 
number, which is already relevant for calculations of the subspace M = — 1. Therefore we have developped 
a drastically improved procedure using cylindrical coordinates which leads to an enormous gain of speed 
such that the computation of the whole Hamiltonian matrix becomes even faster than its diagonalization 
by standard library routines. The derivation of the corresponding powerful formula for the electron-electron 
integral is rather lengthy and complicated: we therefore present only major steps of it in appendices A,B,C 
of the present paper. 



3 Results and discussion 

Throughout the paper we use the notation Vg S+ ^ M^ 1 ^* for a state with spin multiplicity (25+1) and degree 
of excitation u = 1, 2, 3, ... within the subspace of given magnetic quantum number M and z-parity Tl z . The 
index S z will be omitted in obvious cases. The present paper is concerned with the subspaces 1 (— 1) + , 3 (— 1) + , 
x (— 1)~, 3 (— 1)~. The correspondence between our field notation and the common spectroscopic notation 
n|^ +1 Ljvf in field- free space is discussed in ref.|l^] (see table 1 therein). For completeness we mention 
that the 5 1 (±1)~ states correspond to the 6 1 D±\ states whose field free energy is —2. 138982274a. u., and 
analogously the 5 3 (±1)~ states correspond to the 6 3 -D-ti states with a field free energy of —2.013901415a.it. 
(Both values are taken from ref . jig] ) . 
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3.1 Aspects for the selection of basis functions 

For the M = states treated in ref.|L^], we have been able to achieve a considerable accuracy by choosing 
basis sets which can describe the shape of the exact wave function, i.e. include electronic correlation effects. 
The latter become less important with increasing quantum number |M|, and this manifests itself already 
for M = — 1. The reason is that bound two-particle states with nonzero values of M are approximately 
one-particle excitations. Consequently, the electrons are spatially more separated than in a + state, and 
this lowers the correlation energy. Additionally, the cusp problem is also less important for excited states 
M = — 1, and therefore fewer one-particle functions with large values for the nonlinear a and (3 parameters 
are needed. We have exploited these facts and achieved even more accurate results for the M = — 1 states 
than for the M = states. 

In detail the strategy was similar to the M = case. In the case of the (— 1) + subspace we used 260 
optimized one-particle basis functions (for each field strength) for constructing a two-particle basis set of 
dimension N = 3793. The latter number is identical for the singlet and the triplet subspace since it cannot 
occur that a two-particle state contains two identical one-particle contributions which combine to the odd 
total magnetic quantum number M = — 1. This is a principal difference to the case of the + subspace. 
In order to describe angular correlation, we have added also (— 2) + and (— 3) + functions which are paired 
with (+1) + and (+2) + functions, respectively, to build up M = —1. The same scheme was used for the 
contributions of one-particle functions of negative z-parities. There it was sufficient to use the combinations 
(— l) _ /0 _ and (— 2) _ /(+l) _ . In order to describe excitations we added one-particle basis functions with 
quantum numbers m nz = (— 1) + and values k = 1 and ki = 1 (see eqs.(|3|,Q)). The latter have exclusively 
been optimized for a nuclear charge number Z = 1, in contrast to all the other types of basis functions 
which have been optimized for Z = 1 (hydrogen) and for Z = 2(He + ). The reason is again that all bound 
M = — 1 states are one-particle excitations in which the excited electron is associated to the one-particle 
quantum numbers m nz = (— 1) + , and the effect of the nucleus on that outer electron is screened by the inner 
one, thereby giving rise to an effective nuclear charge close to unity. 

For the (— 1)~ subspace we proceeded in a similar manner, the basis set dimension was N = 3671 for 
both singlet and triplet states, built up from a set of 228 optimized one-particle basis functions of type (||). 

In the following, we present and discuss the results for the helium energy calculations. Comparing our 
field- free data with the literature we observe that they are for the subspace (— 1) + /(— 1)~ subspaces more 
accurate than in the case of the + /0~ subspaces. We will as far as possible compare our helium energies 
for finite field strength with the best data available in the literature. 

3.2 Energies for finite field strengths 
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3.2.1 Results for M = — 1 and even z-parity 

a) The singlet states z/ 1 (— 1) + 

For the singlet subspace z/ 1 (— 1) + we present the energies of the ground state and the first four excitations, 
i.e. 1 < v < 5. The energies for the 1) + state are presented in table 1, together with the values given 
in the literature, if available. The accuracy of the field free energy is even higher than the ones of the 1 1 + 
and 1 1 0~ states due to the lower correlation energy for the M = — 1 states. We remark that even though 
the data for the (— 1) + symmetry presented in ref.[|llj have been accurate enough for identifying features 
in the spectrum of the magnetic white dwarf GD229 as absorption edges of helium, we have been able to 
improve the accuracy slightly. 

We observe that for small values of B the total energies of the 1) + state lie slightly lower than for 
B = which is an effect due to the Zeeman energy which is negative for M = — 1. However, for fields 
stronger than B ~ 0.16a.u. the total energies rise drastically with increasing field strength which has its 
origin in the increasing kinetic energy of the electrons which can roughly be estimated by their Landau 
energy amounting to B for both electrons together. 

In order to reveal the internal energetics of the atom we have to subtract such pure and overall field 
effects. For an analysis it is advantageous to subtract even more: one measure for the accuracy of the total 
energies E(B) are their corresponding one-particle ionization energies \E{B) — T{B)\ corresponding to the 
process He — > He + + e~ . The threshold T(B), i.e. the lowest possible total energy for which the system 
He + + e~ can exist possessing the same quantum numbers as the He state in question, is given in the fourth 
column in table 1 (the values for the ionization energies are trivial to compute and are not additionally 
listed in any of the tables). The values T(B) for (— 1) + symmetry are identical to those given in ref.|| for 
the M = symmetry. The reason is that for any z-parity and spin symmetry the lowest M = — 1 state of 
the system He + + e~ is realized by the ionized electron in a Landau state with magnetic quantum number 
m = — 1 and the He + ion in its + ground state: the Landau energy of the electron depends on (m + \m\) and 
is therefore identical for m = — 1 and m = 0. The alternative possibility of associating the value m = — 1 to 
the He + ion possesses a higher energy. 

The energies for the excited states 1) + , 2 < v < 5 are given in table 2. We remark that for finite 
field strengths there are no data about these states available in the literature so far. We observe as expected 
that their total energies rapidly approach the threshold T(B) with increasing excitation. 

The dependence of the one-particle ionization energies on the field strength is shown in figure 1. In 
contrast to the total energies of the 1 1 (— 1) + state its ionization energy depends monotonously on B: the 
outer electron becomes increasingly bound with increasing field strength. This is not in general the case 
for the excited states within the (— 1) + subspace although the ionization energy of the 2 l {— 1) + state is 
monotonous (note the changes in the slope compared to the 1 1 (— 1) + state). The ionization energies of the 
states 3 1 (— 1) + , 4 X (— 1) + and 5 1 (— 1) + exhibit a pattern of several avoided crossings. Those occur in the 
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field strength interval 0.02 < B < 0.2 which is the regime where a rearrangement takes place: for low field 
strengths the states 3 1 (— 1) + and 4 X (— 1) + (field free 4 1 i ? _i and 4 1 P_i) are energetically almost degenerate 
since they both belong to the field free principal quantum number n = 4. This degeneracy is disturbed by 
the magnetic field and completely destroyed for strong fields (see figure 1). 

b) The triplet states z/ 3 (-l)+ 

For the triplet subspace z/ 3 (— 1) + we present the energies for the ground state and the first four excita- 
tions, i.e. 1 < v < 5. The total energies of the states with S z = — 1 are given in table 3 together with data 
existing in the literature for the states 1 3 (-1) + , 2 3 (— 1) + and 3 3 (-l) + . The states 4 3 (— 1) + and 5 3 (-l) + 
have not been present in the literature so far. We observe that our energies are variationally lower than any 
values in the literature apart from only a few exceptions with respect to the energies computed by Jones 
et al [fig] . However, one must take into account that the released-phase quantum Monte Carlo method 
performed by Jones et al leads to statistical error bars (see the corresponding numbers in parentheses in 
table 3) such that none of those energies is significantly stronger bound than our results. This means, either 
we confirm these results or our accuracy is even higher, reflecting the fact that our optimized anisotropic 
Gaussian basis sets excellently describes the wave functions in a magnetic field. Due to the efficiency of our 
method of computation we were able to cover a large number of field strengths. 

Due to the spin shift BS Z the triplet states with S z = — 1 are the most strongly bound ones among the 
three states with S z = 0, ±1. This spin shift causes the 1 3 _ 1 (— 1) + state to cross the low-field ground state 
1 1 + at B ~ 0. 750a.it.. For B ~ 0.750a. -a. the 1 3 _ 1 (— 1) + state is the global ground state of the atom. The 
1 3 _ 1 (— 1) + state is for any field strength lower than the l 3 1 (— 1)~ state (see sec. 3.2.2| ), and it lies also lower 



than the l 3 x + and l 3 1 0~ states for B > 0.750a.n.. The latter ones cross the 1 1 + state at B ~ 1. 112a.it. 
and B ~ 0.994a.?/., respectively. 

Analogously to the singlet case, we show the one-particle ionization energies of the triplet states also 
in figure 1. We observe that the singlet-triplet splitting decreases with increasing excitation. This occurs 
due to the fact that for excitations the spatial separation between the electrons is large, and therefore the 
exchange terms in eqs.([lC|-|i"2"|) are small which causes a small effect of the different signs of the matrix 
elements belonging to singlet and triplet states. 

3.2.2 Results for M = — 1 and odd z-parity 

a) The singlet states u (— 1)~ 

For the singlet subspace u l {—l)~ we present also the ground state and the first four excitations, i.e. 
1 < v < 5. The total energies are given in table 4. Like in the case of the excited 1 (— 1) + states there 
exist no data for finite field strengths in the literature. Our field-free data are in good agreement with the 
literature. In figure 2 we show the one-particle ionization energies analogously to figure 1. The threshold 
T(B) is identical to the case of (— 1) + , + or 0~ symmetry because the energetically lowest way to realize 
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the ionized system He + + e~ with (— 1)~ symmetry is to leave one electron in the + state of a He atom. 
The other electron must then be placed in a Landau orbital with m = — 1 and negative z-parity, which does 
not affect the Landau energy B/2. 

b) The triplet states ^ 3 (-l) 

For the triplet states, several investigations exist in the literature [{15, 18 — pij] , in contrast to the singlet 
case. In table 5 we have listed our total energies and the corresponding values of the literature for the states 
u 3 _ 1 (—l)^ (i.e. S z = —1), 1 < v < 5. Again our results are better than almost all the reference values for 
finite field strengths or are at least comparable to the ones of Jones et al [13] within their statistical error 
bars. For the states 3 3 (— 1) + and 4 3 (— 1) + there exist no data in the literature so far. 

In figure 2 we show the ionization energies. The dashed triplet curves almost coincide with the corre- 
sponding singlet curves. Such a small singlet-triplet splitting is in good agreement with the considerations 



mentioned in sec. 3.2.1 which predict a small singlet-triplet splitting for high excitations. 



4 Transitions 

For the comparison of the energy levels of helium with the spectra of magnetic white dwarfs in general and 
GD229 in particular it is necessary to determine the transition energies from our total energies. Restricting 
ourselves to electric dipole transitions, we have the selection rules AS = 0, AS Z = for the spin degrees 
of freedom in our nonrelativistic approach and AM = 0, AII 2 = ±1 (for linearly polarized radiation) or 
AM = ±1, AII 2 = (for circularly polarized radiation) for the spatial degrees of freedom. Whereas in 



ref.[12] we were already able to present the AM = transitions + <-> CP, we are now able to investigate 
three times as many transitions: firstly the AM = transitions between the (— 1) + and the (— l) - states, 
and additionally two classes of | AM| = 1 transitions between the M = — 1 states and M = states, involving 
positive and negative z-parities, respectively. Altogether our data yield 75 singlet and 70 triplet transitions. 

Due to the fact that the field strengths in the atmospheres of magnetic white dwarfs is not a constant 
but varies by a factor of two for a dipole geometry, transitions which behave monotonically as a function 
of the varying field are smeared out, i.e. are not expected to provide a signature in the observed spectrum. 
However, the transitions whose wavelengths are stationary with respect to the field dependence manifest 
themselves as absorption edges in the observable spectrum if they possess a relevant intensity. We therefore 
give in tables 6 to 11 a complete list of all the stationary points which resulted from our calculated transitions. 

We remark that for a reliable comparison with observational spectra in some cases the finite nuclear 
mass effects have to be taken into account. This requires corrections of our total energies according to the 
scaling relation given in eq.(4) in ref.[^] and in refs-P,[l5]. Due the selection rules AS = 0, AS Z = 0, 
we have for the transition energies the scaling relation AE(M ,^ 2 B) = fiAE(oo,B) - ^/j 2 5AM (here 
Mo = 7344a. u. is the nuclear mass and \x = 0. 999864a. u. is the reduced mass). For AM = transitions 
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the effect is aiways such that the position (field strength) and wavelength of a mass-corrected stationary 
point are related to the corresponding fixed-nucleus result by B(Mq) = [i 2 B{oo) and A(Mq) = ^A(oo). If 
for | AM | = 1 transitions the ratio is small compared to AE it is still possible to do an approximate 
correction for the wavelengths of the stationary points directly with the data presented: We then have 
A(Mo) = ^A(oo)(l + Hjj^\(oo)AM). The stationary points corrected exactly within the scaling relation 
given above can, of course, only be obtained by separately scaling all the values of AE and by interpolating 
them over the grid of scaled field strengths. In the argumentation above we have taken into account the 
normal mass correction terms. The specific mass corrections are expected to be even less significant, in 
particular for stronger fields and excited states. 

Altogether we detected 139 stationary points; several ones among them possess large uncertainties which 
arise mainly due to the interpolation error with respect to the crude grid of field strengths. Most of the 
transitions (— 1) + <-> + , however, are so precise that the corresponding data have already successfully been 
used to explain the absorption edges in the spectrum of the white dwarf GD229|ll]], and these data together 
with the other stationary transitions serve as a good basis for astrophysicists to investigate the spectra of 
unidentified magnetic objects. 

5 Concluding remarks and Outlook 

We have investigated the fixed-nucleus electronic structure of the helium atom in a magnetic field by a fully 
correlated approach. Scaling laws allow to include finite-mass effects. The present work is concerned with 
the energy levels and transitions of the helium states with magnetic quantum number M = — 1. A small part 
of the corresponding data has already been used successfully for identifying the features in the spectrum of 
the white dwarf GD229 with electronic transitions in atomic helium Jll[, which represented one of the goals 
of our work. The enlarged data presented in this paper will now, together with the energy levels provided in 
ref.|l2|], serve as a good starting point for the analysis of observed spectra of magnetic astrophysical objects 
in general. 

The reliability of our wavelength data results from the high accuracy of our energy values which ranges 
between 10 -4 a.u. and 10~ 6 a.u.. This accuracy has become possible due to our appraoch by means of an 
optimized anisotropic Gaussian basis set. Since the spherical invariance is broken by the magnetic field, 
it has been necessary to use Gaussians with different length scales for the longitudinal and transversal 
degrees of freedom. The nonlinear parameters describing these length scales have been determined by the 
requirement to solve optimally the one-particle problem of the H atom or the He + ion in a magnetic field of 
given strength. These optimized one-particle functions have been used to construct configurations in order to 
represent the full fixed-nucleus Hamiltonian and the overlap as matrices separately in each subspace of fixed 
quantum numbers corresponding to the four conserved quantities: the total spin S 2 and its z-component S z , 
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the z-component L z of the electronic angular momentum and the electronic z-parity H z . The corresponding 
generalized eigenvalue problem provided a variational estimation for the energy eigenvalues. 

Atomic energies of helium have been calculated for the ground state and the first four excitations in each 
subspace for M = — 1, i.e. for positive and negative z-parity as well as for singlet and triplet spin symmetry. 
We considered altogether 20 different field strengths < B < WOa.u., i.e. < B < 2.3505 • 10 7 Tesla. 
This series production of data has become by very efficient algorithms for the computation of the matrix 
elements, in particular the electon-electron matrix elements for which we presented an analytical formula 
derived in cylindrical coordinates, thereby making the CPU effort for the calculations for nonzero angular 
momentum affordable: Building up a matrix of dimension about 4000 takes less than one hour CPU time 
on a moderate Silicon Graphics workstation. 

The comparison of energy data with observed spectra of astrophysical objects is possible by searching 
for stationary points of the transitions with respect to the magnetic field strength. The data for M = — 1 
in the present paper yield AM = transitions, and together with the results for M = in ref.^TJ] we have 
also been able to consider AM = 1 transitions. Complete tables of all the detected stationary points for the 
mentioned transitions were given. Energy data for M = —2 and M = —3 are planned to be investigated in 
a future work. 

In order to complete the treatment of bound electronic transitions of helium in a magnetic field, we will 
in the near future also investigate in detail the oscillator strengths of the mentioned transitions as a function 
of the magnetic field. Based on the calculated intensities it will be possible to produce synthetic spectra, 
and a comparison with the observed spectra will give important hints for models of the radiation transport 
and the field configuration in magnetic astrophysical objects. 

Acknowledgements. The Deutsche Studienstiftung (W.B.) and the Deutsche Forschungsgemeinschaft 
(W.B.) are gratefully acknowledged for financial support. 

A Analytical solution to the electron-electron integral 

In the following it is our aim to derive an analytical expression for the electron-electron integral which allows 
its efficient numerical implementation. We emphasize that an efficient treatment of the two-particle integrals 
is essential for the calculations on helium since the number of two-particle matrix elements is N(N + l)/2, 
in contrast to the one-particle matrix elements whose number is only n(n + l)/2 (here N ~ 4000 is the 
dimension of the two-particle Hamiltonian matrix whereas n ~ 200 is the dimension of the underlying one- 
particle basis set, see eqs.(|2],||)). Denoting the two-particle interaction with Vn{^i^2) = ; we have to 
solve the integral 
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where the one-particle orbitals j are of type (||) , obeying the constraints (||J4|) . For the sake of brevity, we 
will in the following use the index notation 7^ := 7i + 7^ for the sum of two indexed quantities. 

The initial step is to apply a Singer transformation [22] in order to remove the Coulomb singularity, 
thereby introducing the new variable u according to | ri ^ r2 | = 77= Jq° du e ~« 2 ( r i- r 2) 2 - Then the integrand of 
this new integration over u decomposes into a transversal and a longitudinal part 

00 

(ij\V n \kl) = A J du I pv {u) • I z (u), (A.2) 


where 

00 00 

I z {u) = J dz x J dz 2 z f>* 3 %'i l e -M-» 2 i*i-'*y , -M (A.3) 
—00 —00 

and 

I pip (u) = J Pl d Pl d Vl p 2 dp 2 dip2 p ^ e -'('».-m t ) We -«, tf j. e -n2(p5^-2 W!mS ( W - W )) ( A4 ) 

x p npjl e ~ i ( m J'~ m i)¥ J 2 e - a iwl 



The longitudinal integral I z (u) is the trivial part of the matrix element ( A.l). For decoupling the particles 1 
and 2 we subsitute Z\ = z\ — b(u)z 2 , z 2 = z 2 where b(u) = p.^ +u i an d d^^j = 1. The exponential factorizes, 
and the power z" zik can be multiplied out, yielding for I z (u) a sum over standard integrals of the Gaussian 
type /f^z^e - ^* 3 w hich can easily be evaluated, giving after a few steps of algebra 

C<^ / x 

I z (u) = 4,g nzijkl J2 y££jA n zik,nzjhPik,Pjr,2C\v) (A.5) 

where the prefactor gi := mod(l,2) reflects the fact that the total z-parity is a conserved quantity. The 
function J is defined by 

_ _ 1 _ 112+1 

J(n u n 2 ,a u a 2 ;v;u) = \a x 1 2 a 2 2 ■ T(x±L) ■ r( "i+"2-u+i ) ■ ( ai a 2 )^ 

nf \ n — Til n 119 — V + l 

x u 2 ^ 1- ^ (1 + iii 2 )V(l + ^n 2 )-^ — (A.6) 

The usual procedure for the treatment of the transversal part I pip would be to represent it in Cartesian 
coordinates since then its decomposition into a sum of products of integrals I x (u) and I y {u) analoguos to 
I z (u) would be technically simple: The exponential would factorize automatically, and the remaining part 
consists of factors like p n p^ m which can be expressed as (x 2 + y 2 ) k (x ± iy)'" 1 ' due to the contraint (|3|) and 
which can trivially be multiplied out. We emphasize that this procedure, however, would yield expressions 
which are very expensive with respect to the CPU time since already for small nonzero maqnetic quantum 
numbers m or excitations k the number of binomic terms gives rise to a large number of integrals over u. 
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In the following we exploit that the described drawback is not at all a property of the integral ( |A,1| ) 
but only of the Cartesian approach. In fact, the many mentioned -u-integrals are not independent from 
each other but the information how to resummarize their results to a more compact expression is lost in 
the lengthy Cartesian algebra. Therefore, it is our strategy to use cylindrical coordinates from the very 
beginning, although the derivation of this condensed analytical formula is technically involved. 

The key for solving the integral I p <p{u) is the following substitution 



Pi = yp? + a?pl - 2apip 2 cos(ip2 ~ fi) P2 = P2 (A.7) 

sin i^i = ^-sin(p 2 - fi) ; cos^i = ^- cos(<p 2 - fi) - a ^ <p 2 = U^ + fl) (A.8) 

Pi Pi Pi 2 

2 

which results from the following ideas (here a(u) = a u +u a )■ Firstly, we use the angular part of I pip {u) for 
exploiting the conservation of L z , and secondly the remaining radial part has to be decoupled similarly like 
in the case of the ^-integration mentioned above. Using in eq.( [A.5 ) the angle <p := ip 2 — ifii and the cyclic 



angle ip := \{<p 2 + <Pi), the integration over ip yields a Kronecker Delta reflecting the conservation of L z : 

oo oo 

I pip (u) = 2tt 6 mi ., mik J Pl d Pl J p 2 d P2 p^ e - a *rt ■ p n 2 P > l e~ a ^ • e "« 2 W+^) 





X J dip e i ( rn i- m k)<P . e -U 2 (-2pip 2 COSip) ^ ^ ( 





Now, the expression (pi cosf) plays a similar role as z\ in eq.( |A.3| ), and therefore we temporarily introduce 
£i = (pi cos f) an d rji = {pi sin ip) as dummy Cartesian variables. For decoupling the particles 1 and 2 
analogously to the treatment of I z (u) we substitute x\ = £i — ap 2 whereas y± = r/i and p 2 = p 2 remain 



unchanged. The variables pi and (p\ given in eqs.( A.7 , A.8 ) are now just the new cylindrical coordinates 



belonging to x\ and y\. In this new representation I pv reads 

2tt oo 







°9 aikajl+a ijkiu 2 

x / dp 2 p 2 P]l e ^+" 2 P2 , (A.IO) 
o 



where the variables p\ and ipi occuring in the factor Pi P%k e^ mi mfc ) vi have to be considered as functions of 
the variables labelled with the symbol ("). We split p\ Plk e *( m *- m fc)</>i = pf ik . p ^+^ e i(mi-m k ) m accor di n g 



to eq.(|3[). This is a powerful step since the relations ( |A.7|JA.8|) can explicitely be solved with respect to those 
two factors. After some algebra we arrive at: 

2fcjfc _ j kik j a 2r 2 +r 3 +r 4 ~2r'i+r 3 +r 4 ~2r 2 +r 3 +r 4 e i(r 3 -r 4 )^ 1 ^ -q^ 



Pi 



.ri r 2 r 3 r 4 . 
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\m,i\+\m k \ i{ mi -m k ) Vl 
Pi e 



\rrii\ 

E 

Mi=0 



\ m i\ 
Pi 



E 



m k\ 1 ^Imil+lmfcl-Mrt ^ifc^l m il+l m fe|-Mifc J{snJ,i-s kt j, k )<pi 



Pk 



P1P2 



(A.12) 



where s, := sgn(mj) are the signs of the magnetic quantum numbers. Inserting the expressions (A.11, A~T2] ) 
into eq.( |A.ld|) , we can perform the (p\ -integration yielding a Kronecker Delta which restricts the summation 
indices to r^ — r^ = —SiHi + s k [i k . The obtained sum is a decomposition of l pip into products of each two radial 
Gaussian integrals completely analogous to the expressions encountered in the evaluation of the integral I z , 



which allows us after several steps to express also I pip in terms of the function J given in eq.( A.6 ): 

\rrii\ 



I ptp (u) 



4"7T 2 5. 



"i hi, ma E 
/j,i=0 



\ m i\ 
Pi 



[Wife I 

E 

Mfe=0 



\m k \ 
Pk 



E 



r 1 +r 2 +r 3 +r 4 = fc ife 
1 "3- r 4 = - s i/ J i+ s *Mfc 



,ri r 2 r 3 r 4 , 



x J ( n Pife + 1 ' n Pji + 1 > a ik-> a jl~, Pik + l + 2ri+r 3 + r 4 ; 



u 



(A.13) 



The final step is to evaluate the it-integration after inserting eqs.( A. 13, A. 5) into eq.(A.2). Considering 

I mj| 



(ij\V n \kl) = 3,27ra6 mU)mij g nMi . M 

fii=0 



\ m i\ 
Pi 



\m k \ 

E 



\m k \ 
Pk 



E 



r 3- r 4 = - s iMi+SA;Mfe 



n r2 n rA 



E 

C=o 



n Z ik 
v2C, 



x J du J(n pik + 1, n pjl + l, a ik , aji; n ik + l + 2r 1 +r 3 + r 4 ; u) ■ J(n zik ,n zjh (3 ik , f3ji;2C;u) (A.14) 


shows that the solution of the electron-electron integral is now reduced to the integration 



K(ni,n 2 ,ai,a 2 ;t>; mi, to 2 , 61, b 2 ; w) 



du J(ni, n 2 , a\, a 2 ; v; u) ■ J (mi, m 2 , b\, 62; iw; u) . (A. 15) 



Up to a constant prefactor, the integrand in eq.( |A.15 ) is of the type g(u) = u 2riu ■ (1 + au 2 ) ra (l + bu 2 ) Tb (l + 
cu 2 ) Vc (l + du 2 ) rd where a = b = jp, c = d = 5^- are fixed real numbers given by the basis functions. 
The exponents depend on the summation indices which enter into the functions J, and n u := ni+m\—v—w is 
running over positive integer values. Further r a := n2 ~ ni and rj, := m2 ~ mi ar e integers whereas r c := v ~ 1 - n ^ 



is integer or half-integer and : = — — 2 mi2 is always half-integer. The substitution a; := 

= id -1 / 2 a;-l/ 2 (l - x )- 3 / 2 cfe leads to 



with 



<?(u) du 



5 / x «»-2(l - x )- f -ab cd -n u -2(l + g aJ ;) r <'(l + %x ) r 'b(l +g c x) rc Cfe (A.16) 



where r a fe c ^ := r a + + r c + and g a := | — 1, := 4 — 1, g c := § — !• I n order to reduce this integral to 
hypergeometric functions it is now necessary to multiply out one of the three last factors. This is possible 
since the basis functions can always be interchanged in such a way that one of the exponents r a , r\, is 
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positive, say r a . Then eq. (3.211) in ref.[23| can be used to obtain 

g(u) du = icT n "-5 jr ( Ta J q s a x B(-r abcd - n u - l n „ + s + 



s=0 



x Fi(n u + s + i,-r 6 ,-r c ,s-r a6cd ;-g 6 ,-g c ) (A. 17) 



where -B ist the beta function and F\ is the Appell hyper geometric function 1 24]. It is defined as double 
series, Fi(a, 6, c; s, y) := £^ =0 fe^y 2*i(a + v, b> , c + 1/; y) where 2 *i(a, 6, c; z) := £^ =0 fegftg 
is the Gaussian hypergeometric function and (a, z/) := ^-fej-^ is the Pochhammer symbol. Including all 
prefactors, we obtain for the integral K: 

K(ni,n 2 ,ai,a2,v;mi,m2,bi,b2,w) 



-ni 



1 _ "12 + 1 



• r(sdhl) • T( n ^-v+i) . ( 0l02 ) 



(A.18) 



1 mio+l 

x 6 1 2 o 2 2 



r(J2±l) • r( m i2-"'+i) • (6i6 2 ) 2 x / y(u) 



where now the integral stands for the expression ( |A.17 ). 

The final result is a real value for (ij \Vn\kl) obtained by inserting eq.( A.l^ ) into eq.( A.14| ): 

\m k \ 



(ij\V n \kl) 



2 - 



mki,mij9n zi j k i ' a n~ 



-m- 



1 "12 + 1 



1 m-to + l l m il /1 

_m i _ 9 a — 9 — \ -> / m 



ji 



E 'T E 



w=0 



Mfe=0 



E 



i-l+r 2 +r3+r 4 = fe ife 
r 3 -r 4 = -s iM j+s fcMft 

^- 2 

- E I 

C=0 

iu gerade 



r\ r 2 r 3 r 4 . 



r(«±i).r(m2^+i).(a W )5 



^ j r(i + c) • r(HHi±i - c) • iPikPnf 



x Fi(n u + s + 1 -r 6 , -r c , s - r abcd ; -q b , -q c ) 



(A.19) 



An overview over the various parameters entering into eq.( [A.19 ) is given in appendix |C] where we briefly 
present the algorithm for the implementation of eq.( A.lSj ). The direct implementation of eq.( A.19| ) yields 
already a factor of 5 to 10 with respect to the increase in speed compared with the common result derived 



in Cartesian coordinates in ref . 1 1 2 ] . However, by systematic exploitation of the symmetries the gain factor 
of CPU can further be raised to between 20 and 40. 



B Symmetry properties of the electron-electron integral 

In the present section we list the different types of symmetries for the expressions encountered in appendix 

E 
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B.l External symmetries 

By an external symmetry we mean a permutation of basis functions which leave Vjj\kl) invariant (up to 
complex conjugation). Such symmetries are firstly (ij\Vn\kl) = (ji\Vn\lk) since the particles are indistin- 
guishable, and secondly we have {ij\Vn\kl) = (kl\ Vn\ij)* due to the hermiticity of Vn, in the case of our 
real matrix elements even (ij\Vn\kl) = (kl\Vn\ij)- 

B.2 Internal symmetries 

The integral {ij\Vn\kl) is also invariant under any manipulation of the basis functions which leave the 



integrals I z {u) and I pv {u) in eq.( A.2 ) invariant separately. This class of internal transformations is larger 
than the class of permutations of basis functions since the parameters entering in I z {u) and I P ip{u) can 
be varied independently. I z (u) is invariant under the simultaneous exchanges i j and k «-> I of indices 
only in the parameters (3 and n z , and I p <p(u) stays unchanged under the analogous interchange of only the 



parameters a, n p and m. Eq.( A.5| ) shows a further internal symmetry of I ptp (u): if the difference rrn — m^ 



remains unchanged the parameters rrii and m/% can be varied arbitrarily as long as the constraint (||) can 
be fulfilled for some k\ and kk, respectively (analogously for m,- and mi). It is an interesting fact that the 
matrix element (ij\Vn\kl) in a subspace with M = mi + mj = + m; possesses the same numerical value 
as some matrix elements involved in subspaces corresponding to some other value of Ml 

The internal symmetries can be exploited for achieving simultaneously n p ^ > n Pik and n Z ji > n zi ^, i.e. 
in the function g(u) (see eqs.( |A.15 , |A.16 )) even both parameters r a and can be assumed to be non-negative. 



A very important consequence of this is that due to the definition of F\ (see paragraph below eq.( A.l?] )) its 



double series reduces to a sum over a finite number of expressions involving the Gaussian hyper geometric 
function 2-^1(^5 b; c; z) which is much simpler to evaluate than F±(a, b, b'; c; x, y). 

B.3 Further symmetries 



Assuming both parameters r a and to be non-negative, eq.( A. 17 ) is also valid after interchanging r a <-> 



(and simultaneously q a <-> qb), representing a new invariance transformation of Vj[\kl). 



C Implementation 



In the following we present the algorithm for the implementation of eq. ( A.19| ) for our calculations on helium 



1. Due the external symmetry (ij\Vn\kl) = (kl\Vn\ij) it is sufficient to compute only the 
upper triangle of the matrix corresponding to Vn. 

2. If rip^ — n Pik < 0, we use the external symmetry {ij\Vn\kl) = (ji\Vjj\lk) and interchange 
globally i <-> j, k «-> /, which achieves r a > 0. 
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3. If then n 



zjl 



n 



zlk 



< 0, we use the internal symmetry of I z (u) and interchange n z 



<-> n 



n zk n zi a s well as /3, <-> /3j, 0k 0h which achieves r b > 0. 
4. We introduce the following abbreviations for the fixed parameters independent on the 
summation indices of eq. ( A.19| ) : 



name 



definition 



nn 
m\ 

r a 

n 
d 

-q a 

-% 



n Pik + 1 
n Pik + n Pji + 2 

n zik 



n zjl n zik 

2 



^ _ PikPjl 



otikiPik+Pil) 



ft, 



1 _ {a-ik+°t]l)PikP]l 
(ftfc+/3ji)"ifc".ji 



type 



integer 
even 

integer 
even 

integer 

integer 
real 
real 
real 
real 



property 



1 < m 

2 < n 12 
< mi 
< m 12 
< r a 

<r b 

1 < d 

-1 < Qa 

-Kq b <0 
-1 < 9c 



5. If then r a > r&, we interchange r a «-> r b , g a «-> g^, thereby minimizing the number of terms 



of the innermost summation in eq.( A.19| ). 

6. Applying the internal symmetry transformation of suitably shifting the values of rrii and 
mk, we replace both rrii and m& by rrii + 2s and + 2s, respectively where s is the 
highest possible integer such that \rrii + 2s | < n Pi and \rrik + 2s| < n Pk . This step is very 
useful for decreasing the number of summations over the indices ri,r2,r^,r^ since raising 
rrii and rrik has the consequence that ki and k k must decrease for keeping n Pi and n Pk 
constant. 

7. We apply the formula ( A.19| ), and the parameters occuring therein depend on the sum- 
mation indices according to the following definitions: 



name 



definition 



type 



property 



v 
w 
n u 
r c 
r,i 

^abcd 



/ij fc + l + 2ri+r 3 + r 4 
2C 

n\ + m\ — v — w 
v— 1— nig 



w—l—mi2 
2 

— «i— mi+ - 1 



integer 

even 
integer 
int. or half-int. 

half-integer 
int. or half-int. 



1< v <n Pik + l 
< w < n zik - g nzik 
9n ztk < n u < n pik + n zik 



1 < Tr. < 



ZJl-l 
2 



_ l+rizik+n Z ji < r , < _ l+9n zik +n Z ji < Q 
2 — d — 2 

-n Pifc - 3 <r abcd < "" p ' fc ~ nz ; fc ~ 3 ~ g " z ^ < 
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8. Although in principle Fi(a, b, b'; c; x, y) is a double power series in the arguments x, y with 
convergence radii |x| < 1 and \y\ < 1 and although in our case even both arguments may 
happen to lie close to 1, we do not have to care about the power series in x. The reason is 
that the corresponding series terminates after 77, terms because the negative integer — 
enters as second parameter argument into F\ in eq.( [A.19 ). The sum terms in F\ each 
involve once the Gaussian hypergeometric function, and since their arguments are related, 
it is possible to use a continued fraction representation for the ratio ^~^^7^g^r~^ i n 
order to establish a recursion law which is stable over the typically < 10 recursions. 
Therefore, for each occurence of F±, it was only necessary to evaluate a single time the 
Gaussian hypergeometric function 2F1, which is not always simple but very efficiently 
possible by means of fast analytical continuation formulas. Several of them are given in 
eqs. (15.3.3-12) in ref.|^5|, a much larger set of such continuation formulas is presented in 
ref.(H. 

We point out that without such a systematic analysis of the electron-electron integral as well as the Appell 
hypergeometric function Fi and in particular the Gaussian hypergeometric function %Fi, the computation of 
many excited helium states for nonzero magnetic quantum number for many different field strengths would 
not have been possible. 
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Captions 



The solid lines correspond to the singlet states, the dashed ones show the triplet states. 



Fig.l. Ionization energies of the singlet and triplet states k l {— 1) + and n 3 (— 1) + , k = 1,...,5, n = 1,...,5. 
The solid lines correspond to the singlet states, the dashed ones show the triplet states. 
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' Fig. 2. Ionization energies of the singlet and triplet states fc (— 1) _ and n 3 (— 1)~, k = 1,...,5, n = 1,...,5. 
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Table 1: Total energies E of the singlet state 1) + , 
one-electron ionization threshold T and, if available, best 
energy values given in the literature, as a function of the 
magnetic field strength B 



B 


E(l 1 (-1) + ) 


literature 


T 


n nnnn 


—9 1 93774 


—9 19384343087 a 


—9 000 

Li .\J\J\J 


n ooos 


—9 1 941 7Q 

— . JL Li ± J. 1 •../ 




— 1 9Q959Q960 


n nod 


—9 195739 

— . -1 .Li ■ I I KjLl 




_1 QQ7qQqn00 


o no8 


— 9 1 97^80 

— . L Li I 00\J 






n 090 

U.UZiU 


— 9 1 ^9^0 




_1 Q8QQ7^001 


0.040 


-2.139008 




-1.979900008 


0.080 


-2.146622 




-1.959600176 


0.160 


-2.148448 




-1.918402804 


0.240 


-2.139976 




-1.876414090 


0.400 


-2.105987 




-1.790105922 


0.500 


-2.077302 


-2.0750 6 


-1.734628064 


0.800 


-1.969560 




-1.561526260 


1.000 


-1.884875 


-1.8800 6 


-1.440989741 


1.600 


-1.590118 




-1.058421519 


2.000 


-1.368986 


-1.3591 6 


-0.788842154 


5.000 


0.619265 




1.456132354 


10.000 


4.500982 




5.609851957 


20.000 


13.010551 




14.47840453 


50.000 


40.339488 




42.45369755 


100.000 


87.671288 




90.43945348 


a Drake et al. (1992) p 


I 



b Larsen (1979) 



Table 2: Total energies E of the excited sing let states i^(-l) + for 2 < v < 5 as a function of 
the magnetic field strength B as well as field free reference values 



B 


E(2*(-l)+) 


E(3*(-l)+) 


£(4'(-l)+) 


£(5 A (-1)+) 


0.0000 


-2.055124 


-2.031253 


-2.031060 


-2.020002 


0.0000 (lit.) 


-2. 055146362" 


-2.031255144° 


-2.03106965° 


-2.020002937° 


0.0008 


-2.055514 


-2.031638 


-2.031422 


-2.020356 


0.004 


-2.056834 


-2.032905 


-2.032089 


-2.021046 


0.008 


-2.057978 


-2.033939 


-2.031432 


-2.020791 


0.020 


-2.058622 


-2.034384 


-2.024218 


-2.017013 


0.040 


-2.053994 


-2.029964 


-2.013028 


-2.003001 


0.080 


-2.038112 


-2.010281 


-1.994622 


-1.982587 


0.160 


-2.004984 


-1.964022 


-1.946526 


-1.937449 


0.240 


-1.970004 


-1.923549 


-1.904394 


-1.893803 


0.400 


-1.893976 


-1.840324 


-1.819308 


-1.808757 


0.500 


-1.843343 


-1.786317 


-1.764444 


-1.753581 


0.800 


-1.680951 


-1.616388 


-1.592658 


-1.581405 


1.000 


-1.565692 


-1.497369 


-1.472743 


-1.461154 


1.600 


-1.194554 


-1.117980 


-1.091462 


-1.079247 


2.000 


-0.930508 


-0.849891 


-0.822478 


-0.809928 


5.000 


1.291512 


1.389207 


1.420194 


1.433857 


10.000 


5.428064 


5.538802 


5.572339 


5.586918 


20.000 


14.279839 


14.403512 


14.439480 


14.455341 


50.000 


42.233602 


42.374094 


42.413014 


42.429203 


100.000 


90.203618 


90.356573 


90.397606 


90.416255 



a Drake et al. (1992) |6j 



3 



Table 3: Total energies E of the triplet states ^ 1 (-1) + (S z = -1) for 1 < v < 5 function of the 
magnetic field strength B. We have also provided the results given in the literature so far. 





£(1 3 (-1) + ) 


£(2 3 (-l)+) 


£(3 3 (-l) + ) 


£(4 3 (-l)+) 


£(5 3 (-l)+) 


D 


this work 


literature 


this work 


literature 


this work 


literature 


this work 


this work 


I"l AAA A 

0.0000 


-2.133149 


-2.133164191° 


-2.058076 


-2.058081084" 


-2.032322 


-2.032324354° 


o no i o cr o 

—2.031253 
—z. 031255108 


o none CT A 

— 2.020550 
— z\0z0551io7 


n nnno. 
0.0008 


-2.134349 


-2.1326 c 


-2.059266 


-2.0588 c 


-2.033487 




O OQO/1 QQ 

— z.Uoz4o8 


— Z.UziOOl 


n nf»/i 
U.UU4 


-2.139109 


— 2.l374 c 


-2. 063816 


—2.0633 


—2.037489 




— z.UoOSoy 


— z.UzoloU 


0.008 


-2.144983 


-2.l437(7) d 


-2.069049 


-2.0689(4) d 


-2.042110 


-2.0407(3) d 


— z. 040 (00 


— Z\Uzoy4U 


U.UzU 


-2.162112 




-2.082183 




-2.054609 




— z.U4osyi 


— z.Uo (oiZ 


n n /i n 
U.040 


-2.189128 


-2.187300 / 


-2.098532 


-2.096752-' 


-2.070637 




— Z.Oool (0 


— Z.U4401 / 


n non 

O.OoO 


-2.238504 


-2.2384(3) 9 


-2.123337 


-2.1231(5) 9 


-2.093337 


-2.0898(3) ff 


— z.O i 48z8 


— z.oozoyo 


n i fin 
O.loO 


-2.325189 




-2.169263 




-2.127206 




O 1 nOQ A Q 

— Z.10yo4o 


— Z.099 /4/ 


n o /i n 
0.240 


— z.4Uzoyo 




— Z.zl4Uo4 


— z.zlUo04 J 


o 1 a. k a no 
— z.l004yz 




o 1 a c c no 
— z. 145502 


1 /i 00 1 
— z.lc>4oZi 


n /i nn 
0.400 


-2.540763 


-2.5366 c 


-2.298465 


-2.2972 71 


-2.241932 


-2.2413' 1 


— z.220053 


— z.z0yz44 


u.ouu 


-2.620021 


-2.6185 6 


-2.348197 




-2.287934 






9 9^4049 

Z .ZidlUIZ 


0.800 


-2.835619 


-2.8356(2)9 


-2.486853 


-2.4869(H) 9 


-2.418154 


-2.4189(2)9 


-2.393403 


-2.381809 


1.000 


-2.965504 


-2.9638 fc 


-2.572178 




-2.499237 




-2.473517 


-2.461576 


1.600 


-3.308774 


-3.3079(7) d 


-2.802301 


-2.8015(8) d 


-2.720071 


-2.6950(12) d 


-2.692304 


-2.679694 


2.000 


-3.508911 


-3.5063 6 


-2.938822 




-2.852078 




-2.823349 


-2.810394 


5.000 


-4.625491 


-4.6173 6 


-3.718569 




-3.613233 




-3.580744 


-3.566622 


10.000 


-5.839475 


-5.8295 e 


-4.582689 




-4.463680 




-4.428595 


-4.413579 


20.000 


-7.440556 


-7.4277 e 


-5.731115 




-5.598915 




-5.561423 


-5.545260 


50.000 


-10.28410 


-I0.2644 e 


-7.777100 




-7.628163 




-7.587802 


-7.571207 


100.00 


-13.10478 


-13.0764 6 


-9.806638 




-9.645518 




-9.603141 


-9.584339 



a Drake et al. (1992) [16] 

b Larsen (1979) [17] 

c Jones et al. (1996) [18] 

d Jones et al. (1997) [19] 

e Ivanov (1994) [20] 

f Thinner et al. (1993) [21] 

9 Jones et al. (1999) [13] 

h ETG Hartree-Fock energies (Jones et al. (1999) [13]) 



Table 4: Total energies E of the states v l {— 1) , 1 < v < 5 as a function of the magnetic field strength B 
well as field free reference values. Reference values for finite field strength are not available in the literature 



far. 



B 


771/1]/ 1 \ — \ 


771/01/ 1 \ -\ 

E(2 L (-1) ) 


771/0I / 1 \ -\ 

E(3 L (-1) ) 


771/ A 1 / 1 \ -\ 

£:(4 i (-l) ) 


£(5 (-1) ) 


0.0000 


-2.055613 


-2.031277 


—2.020014 


-2.020000 


—2.013892 




—2.055620733 


—2.031279846 


—2.020015836 


—2.020000711 


C\ (\ 1 OAAnAA'7/( 

—2.013898227 


0.0008 


-2.056007 


-2.031653 


-2.020375 


—2.020344 


—2.014196 


0.004 


-2.057469 


-2.032717 


-2.021323 


—2.020478 


-2.014483 


0.008 


-2.059049 


-2.033191 


-2.021527 


-2.018934 


—2.013542 


0.020 


—2.062299 


-2.031471 


-2.018372 


—2.012348 


-2.007346 


U.U4U 


O C\RA 9A7 

— 1 


no ^/i /i ^ 

— Z.UZ0440 


— z.uuo 1 iy 


— /.UUU / 00 


— i.yyoz < 


0.080 


-2.060968 


-2.010768 


-1.989904 


-1.979451 


-1.973516 


0.160 


-2.042173 


-1.976616 


-1.951446 


-1.939499 


-1.932936 


0.240 


-2.015792 


-1.939171 


-1.911235 


-1.898364 


-1.891395 


0.400 


-1.951742 


-1.858850 


-1.827207 


-1.813144 


-1.805743 


0.500 


-1.906870 


-1.806041 


-1.772722 


-1.758134 


-1.750501 


0.800 


-1.757739 


-1.638598 


-1.601676 


-1.585991 


-1.577931 


1.000 


-1.649285 


-1.520740 


-1.482091 


-1.465894 


-1.457643 


1.600 


-1.293330 


-1.143718 


-1.101454 


-1.084209 


-1.075561 


2.000 


-1.036760 


-0.876696 


-0.832747 


-0.815025 


-0.806167 


5.000 


1.154910 


1.358659 


1.409032 


1.428522 


1.438083 


10.000 


5.271756 


5.506378 


5.560823 


5.581393 


5.591327 


20.000 


14.109153 


14.370201 


14.427890 


14.449299 


14.459536 


50.000 


42.053685 


42.341109 


42.401827 


42.424005 


42.434562 


100.000 


90.023511 


90.324738 


90.386944 


90.409492 


90.420161 



a Drake et al. (1992) [16] 



Table 5: Total energies E of the triplet states 1) (S z = — 1) for 1 < u < 5 as a function of the magnetic field strength B. We have 

also provided the results given in the literature so far. 





£(i 3 (-in 


E{2\-1)-) 


£(3 3 (-l)-) 






B 


this work 


literature 


this work 


literature 


this work 


literature 


this work 


this work 


0.0000 


-2.055629 


-2.055636309° 


-2.031286 


-2.031288847" 


-2.020020 


-2.020021027° 


-2.020000 
-2.020000711° 


-2.013895 
-2.013901415° 


0.0008 


-2.056824 


-2.056766/ 


-2.032463 


-2.032436/ 


-2.021176 




-2.021148 


-2.014998 


0.004 


-2.061486 


-2. 061428/ 


-2.036727 


-2.036694/ 


-2.025324 




-2.024483 


-2.018484 


0.008 


-2.067062 


-2.067000/ 


-2.041202 


-2.0413(9)°" 


-2.029529 


-2.0266(7)°" 


-2.026940 


-2.021545 


0.020 


-2.082319 




-2.051485 




-2.038378 




-2.032348 


-2.027352 


0.040 


-2.104234 


-2.103728/ 


-2.065465 


-2.057738/ 


-2.048734 




-2.040795 


-2.036278 


0.080 


-2.141017 


-2.1414(3) 9 


-2.090793 


-2.0912(2)9 


-2.069919 


-2.0696(4)9 


-2.059460 


-2.053523 


0.160 


-2.202291 


-2.193460/ 


-2.136663 




-2.111468 




-2.099511 


-2.092944 


0.240 


-2.256006 


-2.238396/ 


-2.179246 




-2.151267 




-2.138381 


-2.131405 


0.400 


-2.352208 


-2.3512' 1 


-2.258990 


-2.2588 h 


-2.227265 


-2.2271 h 


-2.213172 


-2.205759 


0.500 


-2.407521 




-2.306225 




-2.272795 




-2.258169 


-2.250521 


0.800 


-2.559005 


-2.5577(4)9 


-2.438913 


-2.4392(5)9 


-2.401794 


-2.4022(4)9 


-2.386047 


-2.377963 


1.000 


-2.650973 




-2.521134 




-2.482236 




-2.465962 


-2.457680 


1.600 


-2.896192 


-2.8964(7)°' 


-2.744310 


-2.7394(8)°" 


-2.701661 


-2.6436(17)°" 


-2.684304 


-2.675612 


2.000 


-3.040304 




-2.877390 




-2.832986 




-2.815133 


-2.806227 


5.000 


-3.851883 




-3.642446 




-3.591326 




-3.571636 


-3.562005 


10.000 


-4.737490 




-4.494996 




-4.439609 




-4.418795 


-4.408774 


20.000 


-5.902110 




-5.631390 




-5.572603 




-5.550913 


-5.540577 


50.000 


-7.959094 




-7.660647 




-7.598710 




-7.576226 


-7.565565 


100.00 


-9.989376 




-9.677020 




-9.613592 




-9.590738 


-9.579962 



° Drake et al. (1992) [16] 
d Jones ct al. (1997) [19] 
/ Thurner ct al. (1993) [21] 
9 Jones et al. (1999) [13] 

h ETG Hartree-Fock energies (Jones ct al. (1999) [13]) 



Table 6: Overview over all of the stationary points found in the sin- 
glet AM = transitions — 1) + — > 1)~ for fi = 1, 5, ^ = 
1, 5 in the range < i? < 100a. u.. The main part of the errors 
arises due to the interpolation over the relatively crude grid of field 
strengths. For high wavelengths, the finite accuracy of the energy 
values themselves contributes also to error in the energies. 



Component wavelength/A position B/a.u. max/min 

^s+i M (-i) n * 



2 1 


-1)- 


^2' 


-1) + 


15915 ± 50 


0.035 ± 0.003 


min 


2 1 


-1)- 


-+2 1 




16770 ± 85 


0.098 ± 0.007 


max 


2 1 


-1)- 


^3 X 


:-i) + 


13649.7 ± 3.0 


26.06 ± 0.04 


min 


3 1 


-1)- 


-+3 1 




20505 ± 80 


0.0542 ± 0.0004 


min 


3 1 


-1)- 


^3 X 


:-i) + 


37139 ± 20 


0.213 ±0.001 


max 


3 1 


-1)- 


-+4 1 




39258 ± 50 


16.5 ±0.3 


min 


4 1 


-1)- 


-+4 1 




29340 ± 550 


0.08 ±0.01 


min 


4 1 


-1)- 


-+4 1 




76600 ± 600 


0.29 ± 0.01 


max 


5 1 


-1)- 


-+4 1 




25417.4 ± 30 


0.00675 ± 0.00004 


min 


5 1 


-1)- 


-+4 1 




27502 ± 30 


0.03076 ± 0.00001 


max 


5 1 


-1)- 


-+4 1 




18500 ± 1200 


0.059 ± 0.002 


min 


5 1 


-1)- 


-+4 1 




35220 ± 300 


0.23 ± 0.01 


max 


5 1 


-1)- 


-+5 1 




72600 ± 5200 


0.046 ± 0.006 


max 


5 1 


-1)- 


^5 X 




45100 ± 4000 


0.095 ± 0.006 


min 


5 1 


-1)- 


-+5 1 


:-i) + 


208000 ± 12000 


0.290 ±0.015 


max 
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Table 7: Same as Table || but for the triplet AM = transitions 
/i 3 (-l)+ -» ^ 3 (-l)" for /x = 1,...,5, v = 1,...,5 



Component wavelength/A position B/a.u. max/min 

I/ 2S+l M (-l)^ 



2 3 < 


-1)- - 


-3 3 




603000 ± 11000 


0.00320 ± 0.00001 


max 


2 3 


-I) - - 


-3 3 




85000 ± 3200 


0.036 ± 0.002 


min 


2 3 ( 


-1)- - 


-3 3 




139400 ± 700 


31 ±5 


min 


3 3 ( 


-1)- - 


^3 3 




37550 ± 160 


0.0034 ± 0.0009 


max 


3 3 ( 


-1)- - 


^3 3 




17700 ± 300 


0.059 ± 0.001 


min 


3 3 ( 


-1)- - 


-3 3 




32200 ± 60 


0.270 ± 0.006 


max 


3 3 ( 


-1)- - 


^4 3 




39240 ± 200 


0.0045 ± 0.0015 


min 


3 3 ( 


-1)- - 


^ 4 3 




40751 ± 10 


20.0 ±0.5 


min 


4 3 ( 


-1)- - 


^3 3 




12930 ± 240 


0.063 ± 0.003 


min 


4 3, 


-I) - - 


-3 3 




16864 ± 33 


0.215 ±0.004 


max 


4 3 ( 


-1)- - 


^4 3 




32340 ± 270 


0.012 ± 0.0006 


min 


4 3 ( 


-I) - - 


^4 3 




38800 ± 1600 


0.03 ±0.002 


max 


4 3 ( 


-1)- - 


^4 3 




29460 ± 200 


0.081 ± 0.008 


min 


4 3, 


-I) - - 


^4 3 




66600 ± 200 


0.34 ±0.01 


max 


43 ( 


-1)- - 


-5 3 




87175 ± 60 


0.0196 ±0.0002 


min 


5 3 ( 


-l) - - 


^3 3 




11170 ±180 


0.065 ± 0.004 


min 


5 3 ( 


-1)- - 


^3 3 




13470 ± 30 


0.2002 ± 0.0009 


max 


5 3 ( 


-I) - - 


^ 4 3 




23630 ± 40 


0.0098 ± 0.0005 


min 


5 3 ( 


-1)- - 


^4 3 




28200 ± 1000 


0.031 ± 0.003 


max 


5 3 ( 


-I) - - 


^4 3 




21050 ± 520 


0.088 ± 0.01 


min 


5 3 ( 


-I) - - 


^4 3 




32680 ± 340 


0.23 ±0.03 


max 


5 3 ( 


-1)- - 


^5 3 




42200 ± 1700 


0.015 ±0.001 


min 


5 3 ( 


-I) - - 


-5 3 




55790 ± 320 


0.046 ± 0.002 


max 


5 3 ( 


-1)- - 


-5 3 




48400 ± 960 


0.100 ±0.004 


min 


5 3 ( 


-l) - - 


^5 3 




144000 ± 8000 


0.30 ±0.01 


max 


5 3 ( 


-1)- - 


^5 3 




80200 ± 500 


49.4 ± 7 


min 



8 



Table 8: Overview over all of the stationary points found in the 
singlet | AM | = 1 transitions within the subspace of positive z- 
parity, / u 1 (— 1) + — > v 1 + for fi = 1,...,5, v = 1,...,5 in the range 
< B < lOOa.u.. 



Component 

^25+l M (-l) r 






wavelength/ A 


position B/a.u. 


max/min 


P0+ 


-> 




■i; 


+ 


608.1 ±0.3 


0.1672 ±0.0015 


max 


1 1 0+ 




2 X (- 


i; 


+ 


539.95 ± 0.25 


0.0160 ±0.0007 


max 


1*0+ 




3 X (- 


■i; 


+ 


524.83 ± 0.24 


0.017 ±0.002 


max 


1*0+ 




4 X (- 


■i; 


+ 


523.37 ±0.22 


0.004 ± 0.001 


max 


1*0+ 




5 X (- 


■i; 


+ 


516.80 ±0.22 


0.0046 ± 0.0008 


max 


2*0+ 




2 X (- 


-i; 


+ 


5285 ± 3 


0.021 ± 0.003 


max 


2*0+ 




2 X (- 


■i; 


+ 


4812 ± 30 


0.11 ±0.01 


min 


2*0+ 




3 X (- 


■i; 


+ 


4125 ± 3 


0.021 ± 0.003 


max 


2*0+ 




3 X (- 


■i; 


+ 


3417 ± 10 


0.183 ±0.005 


min 


2 1 0+ 




3 X (- 


i; 


+ 


3795.8 ±0.5 


0.802 ± 0.008 


max 


2*0+ 




4 l (- 


■i; 


+ 


4006.2 ± 0.8 


0.004 ± 0.001 


max 


2*0+ 




A l (- 


■i; 


+ 


3015 ± 10 


0.20 ±0.01 


min 


2*0+ 




4\- 


i; 


+ 


3177.7 ± 1.0 


0.661 ± 0.005 


max 


2*0+ 




5 X (- 


■i; 


+ 


3651.7 ±1.4 


0.0050 ± 0.0005 


max 


2*0+ 




5 X (- 


■i; 


+ 


2825 ± 3 


0.22 ±0.01 


min 


2*0+ 






i; 


+ 


2950.4 ± 1.7 


0.637 ± 0.006 


max 


3 1 0+ 






■i; 


+ 


22030 ± 65 


0.034 ± 0.001 


max 


3 1 0+ 






i; 


+ 


12634.5 ± 5.4 


0.177 ±0.002 


min 


3*0+ 






■i; 


+ 


15809 ± 40 


0.00576 ± 0.00008 


max 


3*0+ 






i; 


+ 


8270 ± 150 


0.195 ±0.005 


min 


3 1 0+ 






i; 


+ 


11765 ±6 


1.43 ±0.02 


max 


3*0+ 




5 X (- 


i; 


+ 


11533 ± 26 


0.011 ±0.002 


max 


3*0+ 




5 X (- 


i; 


+ 


7170 ± 20 


0.210 ±0.005 


min 


3*0+ 




5 X (- 


i; 


+ 


9004.4 ± 5.5 


1.08 ±0.02 


max 


4*0+ 




4 X (" 


i; 


+ 


19840 ± 290 


0.004 ± 0.002 


max 


4*0+ 




4 X (- 


i; 


+ 


16190 ± 170 


0.024 ± 0.003 


min 


4*0+ 






i; 


+ 


48900 ± 1400 


0.074 ± 0.009 


max 


4*0+ 




4 X (- 


i; 


+ 


26100 ± 1700 


0.194 ±0.004 


min 


4*0+ 




5 X (- 


i; 


+ 


13283 ± 8 


0.0065 ± 0.0004 


max 


4*0+ 




5 X (- 


i; 


+ 


10490 ± 550 


0.0285 ± 0.0005 


min 


4*0+ 






i; 


+ 


21520 ± 220 


0.092 ± 0.003 


max 


4*0+ 




5 X (- 


i; 


+ 


17100 ± 230 


0.210 ±0.003 


max 
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Table 9: Same as Table |8| but for the triplet \AM\ = 1 transitions 
within the subspace of positive z-parity, ^ 3 (— 1) + — ► z^ 3 + for fj, = 
1,...,5, v = 1,...,5. 



Component 
v 2S+l M {-lf 


z 


wavplpnerth / A 


nosition ff/n a, 


nictx/iiiin 


1°0 + - 


n3^ 

-+ 2°(- 


) 


ZLD^S P, -A- D ZL 


n D9i /i -u n nnm 


IllclX 


1 3n4- 

1°0^ - 


n3/ 

-+ 2 {- 


iH 

/ 


QSOfl -|- 9 

oOoU in z 


n 1 7s -i- n nn9 
u. i / o in u.uuz 


nun 


i 3n+ 

1°0 + - 


o3/ 

-* 3°(- 


- 1 J 




n n9n -i- n nm 
u.uzu in u.uui 


IllclX 


i 3n+ 

1°0^ - 


o3 1 


iH 

1 1 


9C1 C J_ K 


w.zo m u.uz 


mm 


i 3n+ 


o3 / 

-> 3 J (- 


1 1 


Qnon i -u n q 

oUyu. l I u.o 


i 7^7 _i_ n nnzi 
i./o/ in u.uuui 


nicix 


l a + - 


/i3 / 

-> 4°(- 


iH 

1 1 


Q91 n -i- 9 
oziu m z 


U.UUO/D in U.UUUUo 


in&x 


i 3n4- 

l d + - 


/i3/ 

-> 4 a (- 


iH 

/ 


zooo in o 


u.zo in u.ui 


mm 


i 3n4- 

l d + - 


/i3/ 

-> 4 ,3 (- 




9finn zt -i- n ^ 

zuuu.4 m kj.o 


i zird -i- n nno 


m&x 


i 3n+ 


tr3 / 

-> 5 J (- 


■n+ 

1 1 


9Q^Q 9-1-0 7 

zyoy.z in u. t 




nicix 


P0 + - 


r3 / 

-> 5 <5 (- 


iH 
1 1 


99D8 -1- 3 

zzuo m o 


n 9qq -i- n nn^ 
u.zyy in u.uuo 


mm 


i 3r\4- 

1 J + - 


-> 5 d (- 


1 I 


9490 i j. n c 


i Q7Q _i_ n nno 
i.o / y in u.uuy 


in3iX 


n3r\+ 

2 d + - 


-> 3 J (- 


1 1 


1 7RQD -1- 8D 
i < oyu in ou 


n nzL7 -i- n nn9 
u.uqi in u.uuz 


nicix 


2 3 + - 


o3 / 

-> 3 J (- 


) 


1 1 n^s 9 -i- 9 n 
iiuoo.z in z.u 


n 938 -i- n nn^ 
u.zoo in u.uuo 


mm 


n3r\+ 


/i3/ 

-> 4< 5 (- 


iH 
1 1 


1 9QQD -1- 39 
izyyu m oz 


u.uuoo in u.uuu^t 


nicix 


n3r\+ 

2 J + - 


-> 4 a (- 


1 1 


741 n -1- zin 
I ^iu m ^±u 


n 97 + n (19 
u.z i in u.uz 


mm 


r»3r*4- 


^ 3 / 

-> 4 ,3 (- 


iH 

/ 


±u^±oo in y 


9 QQ _|_ n D7 

z.yo in u.u t 


nicix 


2 3 0+ - 


-5 3 (- 


iH 
J 


Q71 1-1-13 
y ( 1 1 in lo 


n ni 3 -i- n nn3 
u.uio in u.uuo 


max 


2 3 0+ - 


-5 3 (- 


iH 

J 


ozou in 4u 


n 9Q + n (ii 
u.zy in u.ui 


mm 


2 3 0+ - 


-5 3 (- 


i) + 


8031.8 ± 1.7 


2.14 ±0.08 


max 


3 3 0+- 


^4 3 (- 


i) + 


20450 ± 120 


0.007 ± 0.001 


max 


3 3 0+- 


+ 4 3 (- 


i) + 


14850 ± 300 


0.035 ± 0.003 


min 


3 3 0+- 


-4 3 (- 


i) + 


37220 ± 680 


0.094 ± 0.003 


max 


3 3 0+- 


+ 4 3 (- 


i) + 


26080 ± 250 


0.281 ± 0.009 


min 


3 3 0+- 


^5 3 (- 


i) + 


13320 ± 3 


0.0078 ± 0.0009 


max 


3 3 0+- 


^5 3 (- 


i) + 


11150 ±560 


0.032 ± 0.003 


min 


3 3 0+- 


-5 3 (- 


i) + 


20300 ± 1000 


0.114 ±0.007 


max 


3 3 0+- 


-5 3 (- 


i) + 


16030 ± 210 


0.29 ±0.01 


min 


3 3 0+- 


^5 3 (- 


i) + 


24864 ± 11 


4.19 ±0.08 


max 


4 3 0+ - 


^4 3 (- 


i) + 


296800 ± 8500 


0.009 ± 0.001 


max 


4 3 0+ - 


^4 3 (- 


i) + 


297000 ± 8500 


0.26 ±0.02 


max 


4 3 0+ - 


-5 3 (- 


i) + 


42500 ± 1300 


0.017 ±0.002 


max 


4 3 0+ - 


-5 3 (- 


i) + 


82600 ± 3300 


0.15 ±0.02 


max 


4 3 0+ - 


^5 3 (- 


i) + 


46400 ± 2000 


0.29 ±0.02 


min 


4 3 0+ - 


^5 3 (- 


i) + 


240400 ± 1400 


8.3 ±0.5 


max 
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Table 10: Overview over all of the stationary points found in the 
singlet | AM | = 1 transitions within the subspace of negative z- 
parity, fi (— 1) _ — > z^O - for /j, = 1, 5, v = 1, 5. 

Component wavelength/A position B/a.u. max/min 

i,2S+l M (-l) n * 



2 1 




+ i*o- 


5056 ± 22 





0088 ± 0.0002 


max 


3 1 


-1)- - 


- ^o- 


4477 ± 18 





0070 ± 0.0003 


max 


3 1 


-1)- - 


-4 2 1 _ 


13875 ± 42 





0105 ± 0.0001 


max 


3 1 


-1)-- 


+ 3 X 0- 


56300 ± 1500 





0139 ± 0.0015 


max 


3 1 


-1)-- 


-» 3 X 0~ 


50841 ± 100 





0387 ± 0.0035 


min 


4 1 


-1)-- 


+ 1*0- 


4442 ± 20 





0022 ± 0.0006 


max 


4 1 


-1)- " 


+ 2*0" 


13327 ± 75 





0025 ± 0.0007 


max 


4 1 


-1)- - 


-» 3^- 


44050 ± 570 





0027 ± 0.0009 


max 


5 1 


-1)- - 


- ^o- 


4188 ± 16 





0030 ± 0.0008 


max 


5 1 


-1)-- 


+ 2 1 _ 


11292 ±32 





0035 ± 0.001 


max 


5 1 


-1)-- 


-» 3*0" 


27640 ± 30 





004 ± 0.001 


max 


5 1 


- 


+ 4 1 0~ 


48860 ± 280 





195 ± 0.003 


max 


5 1 


-i)- - 


-» 4 X 0~ 


45700 ± 300 


5 


2 ±0.2 


min 


5 1 


-i)- - 


-» 5 1 0" 


98980 ± 600 





0065 ± 0.0003 


max 


5 1 


-i)- - 


5*0- 


83470 ± 420 





0221 ± 0.0009 


min 



Table 11: Same as 
the subspace of neg 
H = 1,...,5, v = l, 



Table |Ig but for the triplet \AM\ 
;ative z-parity, transitions /U 3 (— 1)~ 
...,5. 



1 within 
i/ 3 0~ for 



Component 



wavelength/A position B/a.u. max/min 





.l)n« 








1 3 ( 


-1)-- 


-> 1 3 0~ 


7143 ± 30 


0.1103 ±0.0007 


max 


2 3 ( 


-1)-- 


-» 1 3 0~ 


4574 ± 15 


0.0087 ± 0.0003 


max 


3 3 ( 


-1)-- 


-» 1 3 0~ 


4095 ± 10 


0.0070 ± 0.0003 


max 


3 3 ( 


-1)-- 


-» 2 3 0~ 


12680 ± 22 


0.00965 ± 0.00006 


max 


3 3 ( 


-1)- - 


-» 3 3 0" 


51500 ± 1300 


0.013 ±0.001 


max 


3 3 ( 


-1)-- 


-» 3 3 0~ 


43530 ± 400 


0.046 ± 0.004 


min 


4 3 ( 


-1)-- 


-» 1 3 0~ 


4063 ± 12 


0.0024 ± 0.0003 


max 


4 3 ( 


-1)-- 


-» 2 3 0~ 


12263 ± 25 


0.00224 ± 0.00007 


max 


4 3 ( 


-1)-- 


-» 3 3 0~ 


40265 ±110 


0.00264 ± 0.00001 


max 


4 3 ( 


-1)-- 


-» 3 3 0~ 


23021 ± 40 


0.079 ± 0.007 


min 


4 3 ( 


-1)-- 


-» 3 3 0~ 


23532 ± 50 


0.17 ±0.01 


max 


5 3 ( 


-1)-- 


-» 1 3 0~ 


3850 ± 9 


0.0033 ± 0.0009 


max 


5 3 ( 


-1)-- 


-» 2 3 0" 


10505 ± 20 


0.0035 ± 0.0007 


max 


5 3 ( 


-1)- - 


-» 3 3 0~ 


26225 ± 70 


0.0050 ± 0.0005 


max 


5 3 ( 


-1)-- 


-» 4 3 0~ 


44662 ± 30 


0.24 ±0.02 


max 


5 3 ( 


-1)-- 


-» 5 3 0~ 


91500 ± 1200 


0.0062 ± 0.0003 


max 


5 3 ( 


-1)-- 


-» 5 3 0~ 


2950 ± 530 


0.022 ± 0.001 


min 
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